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We construct and discuss a semi-rational, multi-parametric vector solution of coupled nonlinear 
Schrodinger equations (Manakov system). This family of solutions includes known vector Peregrine 
solutions, bright-dark-rogue solutions, and novel vector unusual freak waves. The vector freak (or 
rogue) waves could be of great interest in a variety of complex systems, from optics to Bose-Einstein 
condensates and finance. 
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Introduction. Extreme wave events, also referred to as 
freak or rogue waves, are mostly known as an oceanic 
phenomenon responsible for a large number of maritime 
disasters. These waves, which have height and steep- 
ness much greater than expected from the sea average 
state pj , have recently become a topic of intense research. 
Freak waves appear both in deep ocean and in shallow 
waters Q . In contrast to tsunamis and storms associated 
with typhoons that can be predicted hours (sometimes 
days) in advance, the particular danger of oceanic rogue 
waves is that they suddenly appear from nowhere only 
seconds before they hit a ship. The grim reality, however, 
is that although the existence of freak waves has now been 
confirmed by multiple observations, uncertainty remains 
on their fundamental origins. This hinders systematic 
approaches to study their characteristics, including the 
predictability of their appearance 0] . 

In fact, research on rogue waves is in an emerging state 
[E B EI • These waves not only appear in oceans Q but 
also in the atmosphere Q, in optics [1, 0, in plas mas 
[Til ] . in superfluids, in Bose-Einstein condensates jl2| and 
also as capillary waves The common features and 

differences among freak wave manifestations in their dif- 
ferent contexts is a subject of intense discussion Q. New 
studies of rogue waves in any of these disciplines enrich 
their concept and lead to progress towards a comprehen- 
sive understanding of a phenomenon which still remains 
largely unexplored. A formal mathematical description 
of a ro gue wave is provided by the so-called Peregrine 
soliton [14j. This solitary wave is a solution of the one- 
dimensional nonlinear Schrodinger equation (NLSE) with 
the property of being localized in both the transverse and 
the longitudinal coordinate: thus it describes a unique 
wave event. This solution is also unique in a mathemat- 
ical sense, as it is written in terms of rational functions 
of coordinates, in contrast to most of other known so- 
lutions of the NLSE. Recent experiments have provided 
a path to generating Peregrine solitons in optical fibers 
with standard telecommunication equipment [Tol ]. The 
further experimental observation of Peregrine solitons in 
a water tank [16[ indicates that they can also describe 
rogue waves in oceans. The Peregrine soliton is not the 
only fully localized waveform [l7| . In fact, there is an in- 



finite hierarchy of rational solutions of the NLSE which 
enjoy the same property [H|-[2l[. 

In a variety of complex systems such as Bose-Einstein 
condensates [HI, optical fibers 23 1, and financial sys- 
tems [HI HH , several variables rather than a single wave 
amplitude need to be considered. For instance, in the fi- 
nancial world it is necessary to couple cash to the value of 
other assets such as shares, bonds, options, etc., as well 
as to consider all correlations between these variables. 
The resulting systems of equations may thus describe ex- 
treme waves with higher accuracy than the single NLSE 
model. Approaches to rogue wave phenomena involving 
multiple coupled waves are the coupled Gross-Pitacvskii 
(GP) equations [26j | and the Manakov system (or vector 
NLSE) [23]. Indeed, vector rogue waves of GP equations 
and the Manakov system have been recently presented 

In this Letter, we construct and discuss a novel semi- 
rational, multi-parametric vector solution of the Man- 
akov system. For special parameter values this solution 
reproduces known vector rogue waves (such as the vec- 
tor Peregrine soliton [22l and bright-dark-rogue waves 
Our treatment below goes as follows. We give 
the essential Darboux dressing transformation to con- 
struct freak solutions of the Manakov system. We present 
the multi-parametric, semi-rational deterministic freak 
waves. Moreover, we discuss their experimental feasi- 
bility in nonlinear optics. 

Darboux dressing technique. Waves are assumed to be 
modeled by the dimensionless vector coupled nonlinear 
Schrodinger equations (VNLSE) or Manakov system: 
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where u^\x,t), u^ 2 \x,t) represent the wave envelopes 
and t, x are the longitudinal and transverse coordi- 
nates, respectively. Each subscripted variable in Eqs. 
([T|) stands for partial differentiation. It should be 
pointed out that the meaning of the dependent variables 
vS^'ix, t),u^(x, t), and of the coordinates t,x depends 
on the particular applicative context (f.i. plasma physics, 
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nonlinear optics, finance). Note also that Eqs. (JT|) re- 
fer to the self-focusing (or anomalous dispersion) regime. 
Eqs. ((l) are intcgrable: the associated Lax pair is 

$> x = (ika + Q)# , * t = [2ik 2 a + 2kQ + ia(Q 2 - Q x )]%, 

.(2) 

where = k) is a 3 x 3 matrix solution, 

k is the complex spectral variable, the matrix a = 
diagjl , —1,-1} is constant and diagonal, and Q = 
Q(x,t) is the 3x3 matrix 
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Formula ((6j) shows that, if the matrix A (and therefore f2) 
possesses three linearly independent eigenvectors, then 
the vector Z(x,t) is a linear combination of exponential 
functions of x and t. Therefore the solution ([5]) cannot 
be rational or semi-rational. Since only the exponential 
function of a nilpotent matrix is polynomial, one has to 
find those particular values of /3 (see ([7])) such that the 
matrix A (and therefore f2) is similar to a Jordan matrix. 
Indeed this happens if and only if p = ±ia. By taking 
f.i. ft — ia, in this way we arrive to the following semi- 
rational solution of the VNLSE Eqs. (P) 
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The starting point here is the construction of the solution 
representing one soliton wave nonlincarly superimposed 
to the following plane wave background solution of Eqs. 
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where a\ and 02 are arbitrary parameters which, with 
no loss of generality, are taken real. Moreover the fre- 
quency u! reads as oj — a 2 where, from now on, we set 
a = \/ a\ + a 2 . The Darboux method to obtain such 
one-soliton solution (x, t), (x, t) is well known, 
therefore we omit detailed computations, limiting our- 
selves to list the essential few steps. Our results have 
been obtained by following t he g eneral formulation and 
construction as presented in [29[ (the interested reader 
may find additional references quoted there). 

The chain of calculations ends up with the following 
general formula [29| 
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The constant parameter j3 is complex (with non 
vanishing imaginary part), while the functions 
((x, t) , z^'(x, t) , z^ 2 )(x, t) are the components of a 
generic 3-dimensional vector solution Z(x,t) of the Lax 
pair of equations ([2"]). which corresponds to the spectral 
parameter k — f3 and to the background solution (0}. 
Thus, if Zq is an arbitrary complex 3-dimcnsional 
vector, Z(x, t) reads as 
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cxp(iAx — iSlt)Zo, 
(6) 



with the following notation: L = | — 8u} 2 t 2 — 2a 2 x 2 + 
8iwt + \ f\ 2 e 2ax , M = if (ax - 2iwt - I) e ( aa: + 8 "*) , and 
B = ^+8uj 2 t 2 + 2a 2 x 2 + \f\ 2 e 2ax , where / is a complex ar- 
bitrary constant. It should be remarked that the dressing 
construction of the vector rogue wave ([8]) has introduced 
as arbitrary parameters the three complex components of 
the vector Zq, see (O and (J6j> . However only the complex 
parameter / is left essential out of these components, 
since the other parameters enter as trivial translations 
of the coordinates x , t, which have been set to zero for 
simplicity. The two other real parameters ai,a 2 origi- 
nate instead from the naked solution, namely from the 
background plane wave ((4]) . We note also that the depen- 
dence of L, M and B (see jS])) on x, t is both polynomial 
and exponential only through the dimensionless variables 
ax and cot = a 2 t. Moreover the vector solution © turns 
out to be a combination of the two constant orthogonal 
vectors (ai , a.2) T and (02 , — ai) T . 

Vector semi-rational rogue waves. Setting / = im- 
plies M = 0: in this particular case the expression © 
yields the trivial vector generalization of the rational 
Peregrine solution [THHH. In this case u^(x, t) is merely 
proportional to u^(x,t). For future reference we note 
that the amplitude \u^\x,t)\ is picked at x = with the 
maximum value 3|aj| at t = 0. 




where A and Q are the constant matrices 



FIG. 1: Deterministic vector freak waves envelope distribu- 
tions Ivr 1 ' (x, i)\ and |u' 2 '(x,t)| of @. Here, ai = l,a2 = 

0.1, / = 0. 



A = 



j3 ia\ ia2 
-ia\ —(3 
-ia2 —f3 



n = 



-A 2 -2/3A + /3 2 + 2a 2 . 

(7) 



If instead / 7^ 0, the Peregrine bump coexists and in- 
teracts with a pulse which propagates with non constant 
speed, and, depending on the value of |/|, may have dif- 
ferent looks. In order to better describe this behavior wc 



note that the ratios L(x,t)/ B(x,t) and M(x,t)/ B(x,t) 
which appear in ((5J) describe asymptotically as t — > ±00 
a dark and a bright pulse, respectively. Thus each wave 
component «»' (x, t) is generically a mixture of a dark 
and a bright pulse. Leaving aside the detailed analytic 
description of the solution (J8} at intermediate times, we 
limit our present analysis to the large time behavior. The 
pulse motion, for each individual dark and bright contri- 
bution, asymptotically reads as 



x = £(t) -> x + -ln(w|f|) + 



ln 2 (|t|) 
t 2 



±00, 



(9) 

where xq = (1/a) ln(2 v / 2/|/|). This implies that this 
pulse goes to infinity where it "stops" since its velocity 
slowly vanishes as d£ / dt 1 / (at) . The shape of the dark 
and bright contributions at large times as a function of 
the parameter \ which measures the displacement from 
the peak position, takes the expected form as t — > ±00 



L w x Af r-( ft\ e iut 

-^tanhW,-^-^^-^. (10) 

The superposition of the dark and bright contributions 
in each of the two wave components \u^\ may give rise 
to complicated breather-like pulses. These results are 
well represented in Figs. [2]|4j The single contributions of 
the dark shape L/B and bright shape M/B are better 
displayed when f.i. ai = 0. In this case typical distri- 
butions \u^'(x,t)\, \u^(x,t)\ are displayed in Figs. [2j 
|3l Figure [2] shows a vector dark-bright soliton together 
with a single Peregrine soliton. Decreasing the value of 
|/|, Peregrine and dark-bright solitons separate. By in- 
creasing I /|, Peregrine and dark-bright solitons merge 
and the Peregrine bump cannot be identified while the 
resulting dark-bright pulse apprears as a boomeron-type 
soliton (see Fig. [3]), i.e. a soliton solution with a time- 
dependent velocity [3(1 HJ ■ Note that the solution flSJ 
includes as a special case the bright-dark-rogue wave so- 
lution that was reported in (28|. 




FIG. 2: As in Fig. [Q with / = 0.1, ai = 1, a 2 = 0. 

Finally our formula ([5J), if all parameters /, 01,02 are 
non vanishing, describes the dynamics of a breather- 
like wave resulting from the interference between 
the dark and bright contributions. Distributions 
(x , t)\ , (x , t)\ which are typical of this general 



FIG. 3: As in Fig. [T] with / = 10, 01 = 1, a 2 = 0. 



case are displayed in Fig. |4j Again, by decreasing |/| 
Peregrine and breather solitons separate while Peregrine 
and breather solitons merge, with boomeronic behavior, 
if instead |/| increases. 





FIG. 4: As in Fig. [TJ with / = O.li, 01 = 1.2, a 2 = 1.2. 

These results provide evidence of an attractive inter- 
action between the dark-bright wave and the Peregrine 
rogue wave. The observed behavior may also be inter- 
preted as a mechanism of generation of one rogue wave 
out of a slowly moving boomeronic soliton. 

Let us discuss the experimental conditions for the ob- 
servation of the vector, semi-rational freak solitons. Non- 
linear optics is a fertile ground to develop the knowl- 
edge of the phenomenon of vector freak or rogue waves. 
As first scenery, consider the propagation of arbitrar- 
ily polarized optical pulses in a weakly dispersive and 
nonlinear dielectric. In fact, Eqs. (|TJ) apply to a Kerr 
medium with the clcctrostrictive mechanism of nonlin- 



earity 3^ , as well as to randomly birefringent fiber optic 
transmission links [H, [3(| . Indeed, the use of the po- 
larization degree of freedom for doubling the capacity of 
long-distance fiber-based transmission systems has been 
currently widely adopted by means of the technique of 
polarization multiplexing. To be specific, we consider 
the transmission at the 40 Gbit/s rate of a train of dark- 
bright solitons, dark in one polarization, and bright in the 
orthogonal polarization. Fig. [5] shows that a Peregrine 
soliton is generated at 800 km, and it attracts a dark- 
bright soliton. In this example, we numerically integrated 
Eqs. (U} for properly rescaled wave envelope amplitudes 
Ey, Ex and rescaled coordinates, with initial conditions 
two dark-bright solitons plus a small noise seed. We used 
a fiber nonlinear coefficient of 1.3 km _1 W _1 , the anoma- 
lous average fiber dispersion of 0.1 ps 2 /km, and a dark- 
bright full width at half maximum of 8.25 ps (33% of the 
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FIG. 5: Numerical transmission of two 50 ps spaced dark- 
bright solitons in optical fibers, y-polarized dark waves (Ey), 
and x-polarized bright envelopes (Ex)- 

25 ps bit period); the peak power of the two polariza- 
tions is equal to 3 mW and 6 mW, respectively. As sec- 



ond scenery, we may consider incoherently coupled pho- 
torefractive spatial waves in strontium barium niobate 
(SBN). Modulation instability and the existence of un- 
stable dark-bright pairs (first steps in demonstrating vec- 
tor Peregrine waves and dark-bright- Peregrine dynamics) 
have been already demonstrated in SBN [32[. Set-ups 
proposed in Ref. [32l . [33| can be exploited to observe and 
characterize spatial vector rogue waves in SBN. 

Conclusions. Here we have analytically constructed 
and discussed, a multi-parametric vector freak solution 
of the vector NLSE. This family of exact solutions in- 
cludes known vector Peregrine (rational) solutions, and 
novel freak solutions which feature both exponential and 
rational dependence on coordinates. Because of the uni- 
versality of the vector NLSE model ([TJ, our solutions 
contribute to a better control and understanding of rogue 
wave phenomena in a variety of complex dynamics, rang- 
ing from optical communications to Bosc-Einstein con- 
densates and financial systems. 
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